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Abstract 

Gravitational perturbations in the geometry of asymptotically flat black hole 
solution in the deformed Hofava-Lifshitz(HL) gravity are considered and the 
associated quasinormal mode(QNM) frequencies are evaluated using WKB 
method. Black holes are found to be stable against small perturbations. We 
also find a significant deviation in the behavior of QNMs in HL theory for low 
values of HL parameter u, from that in the standard General Relativity(GR). 
QNMs are long lived in Hofava theory compared with the usual Schwarzschild 
black hole case. 
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1. Introduction 

One of the greatest goals in theoretical physics is to realize a quantum 
version of gravity. Attempts are being made in the past to explain gravity 
in the framework of quantum field theory but none are satisfactory. Re- 
cently a renormalizable theory of gravity in 3+1 dimensions was proposed 
by Hofava[l[, inspired from the Lifshitz theory in solid state physics, now 
known as Hofava-Lifshits theory. The theory is a potential candidate of 
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quantum field theory of gravity. It assumes a Lifshitz-like anisotropic scal- 
ing between space and time at short distances, characterized by a dynamical 
critical exponent z = 3 and thus breaking the Lorentz invariance. While in 
the IR limit it flows to z = 1, retrieving the Einstein's GR. Even though 
the discussions on some issues with the theory are going on0, 0, 0, 0], as 
a new theory, it is interesting to investigate its various aspects in parallel. 
Cosmological implications 0,0, 0, 0, Eli 11 1 and astrophysically viable tests 
of the theory [12J, [13j, [14|, [15| had been studied. Since the theory has the 
same Newtonian and post Newtonian corrections as those of GR, systems 
of strong gravity, like black holes, are needed to get observable deviation 



from the standard GR. The black hole solutions in HL theory were stud- 



and some of its properties were investigated 



Stability of black holes in HL theory is the topic of this paper. There are 
strong indirect evidences for the existence of black holes, and if the black holes 
are found to be unstable in the new theory it may question the reliability of 
the theory itself. 

Gravitational perturbation and stability studies of black holes were orig- 
inally formulated by Regge and Wheeler 29, 3(|. Considering small pertur- 
bations in the space time outside the horizon and assuming harmonic time 
dependence to this perturbation, seek the solutions of the perturbation equa- 
tion by imposing the boundary conditions that the perturbations will vanish 
at the horizon and at spatial infinity. If it admits solutions with damping 
amplitudes (frequencies with negative imaginary part) then the black hole is 
stable, otherwise the perturbations will grow in time and the gravitational 
waves will carry away energy of the system leading to an unstable black hole. 
This complex frequency with negative imaginary part is called quasinormal 
modes and its existence proves the stability of the black hole. Also features 
of the black holes in the new theory would be imprinted on the QNMs and 
if the gravitational detectors can observe them it can test the validity of the 
theory. 

The paper is organized as follows. In the following section we briefly 
review asymptotically flat black hole solution in modified HL theory and 
derive the radial equation for perturbation in that black holes spacetime. 
In sectioii|2] we evaluate the QNMs. The main conclusions are discussed in 
section^ 
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2. Gravitational perturbation around black hole in Hofava theory 



The IR vacuum of pure HL gravity is found to be anti-de Sitter [lfj, [18 . 
Even though HL gravity could recover GR in IR at the action level for a 
particular value of the parameter A = 1, there found a significant difference 
between these black hole solutions and the usual Schwarzschild AdS. The 
asymptotic fall-off of the metric for these black hole solutions is much slower 
than that of usual Schwarzschild AdS black holes in GR. Meanwhile Keha- 
gias and Sfetsos[17j could find a black hole solution in asymptotically flat 
Minkowski spacetimes applying deformation in HL theory by adding a term 
proportional to the Ricci scalar of three-geometry, fx^R^ and then taking 
the limit Kw ~~ ^ 0- This will not alter the UV properties of the theory but it 
does the IR ones leading to Minkowski vacuum analogous to Schwarzschild 



spacetime in GR. The deformed action in this particular limit is given as [17 



S = f dtd^ 9 N^(K l3 K- - \K 2 ) - ^QjC* + g^4 3) vX 3) ' 

_ k ^ R m R m + 8( f_^ A) 1 V 3) ) 2 + ^ (3) }, (i) 



with second fundamental form defined as: 



Kij = ^ {jto ~ VtNj - VjM) , (2) 



and Cotton tensor C*- 7 , 



C ij = e lkl V k - ^R {3) 5^ . (3) 

Considering the metric ansatz for a static, spherically symmetric solution, 
ds 2 = -N(r) 2 dt 2 + fir^dr 2 + r 2 dil 2 , (4) 
the Lagrangian after angular integration reduces to 



C 



fcV N 
8(1-3A)77 



(2A - 1)^1L _ 2 \L-±f + ±-±f> 2 -2u(l-f- rf) 

7* IT Zi 



(5) 

In the IR limit A = 1 of the theory, one gets asymptotic flat solution 
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N(r) 2 = f(r) = 1 + ur 2 - y/r(u; 2 r 3 + AuM), (6) 
with horizons at, 

r± = M ( 1± \/ 1 -2^)- (7) 

where M is the integration constant related to ADM mass. In order that 
the singularity at r = not to be naked it requires uM 2 > ~. In the limit 
ujM 2 ^> 1 this solution recovers Schwarzschild solution in conventional Gen- 
eral Relativity. Metric with coefficients given in Eq.([6]) is our initial 
time-independent equilibrium system up on which we consider small pertur- 
bation with harmonic time dependence e~ tkt . Where k is the frequency. 
Now at late time, we want to see whether this frequency has an imaginary 
part which makes the perturbations to grow exponentially in time(unstable) 
or leads to its damping(stable). 

Equation governing the perturbation is the Einstein's equation for this 
perturbed spacetime, 

R P v{g + h) = 0, (8) 

where R^^g + h) is the Ricci tensor computed from the total metric g^v+h^. 
Since the perturbation is small, taking terms up to linear in h^, above 
equation can be expanded as, 

R^ig) + 8R pv {h) = 0, (9) 

where R liv {g) is the Ricci tensor computed from the unperturbed metric g^ v 
which will vanish and the equation becomes, 

SR^{h) = 0. (10) 



Variation of Ricci tensor can be found from the expression 31 



5R„ = -5T^ + 5T% ;u , (11) 
where the variation of affine connections is 

^r^ 7 = -jg av {hp v -^ + hry U .p - hp r>v ), (12) 

Any arbitrary perturbations can be decomposed in to normal modes, since 
the background we are considering is spherically symmetric. For any given 
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value of the angular momentum L, associated with these modes, there are 
two classes of perturbations. Even, (— 1) L and odd, (— 1) L+1 parity pertur- 
bations. Here we consider odd parity case. The canonical form of odd wave 
perturbations in Regge- Wheeler gauge [291] is 



U [LV 



h (r) 

hi{r) 



ho(r) ht{r) 



d 



e- M (sin(e)-P L (cos(e))). (13) 



Substituting Eq. (jl3jl in Eq. (lT0l) and separating the angular and radial 
parts we get the following radial equations: 



ikhn(r) d r „, . „ _„ 

~~f + ~dr [/ l(r)] = ' 23 ' 



(14) 



ik ( dho 



2h r 



f(r) \ dr r 



2u(r 



-y/ra;(4M+r 3 aj) ' 



/ 



0, from 6R 13 . (15) 



Defining dr* = jdr and Q = ^y-, the above two equations can be com- 
bined to get a second order equation by eliminating ho(r): 

d 2 Q 



where 



eff 



f 



drl 
L(L + 1) 



+ (k 2 - V eff )Q = 0. 



6a; 1 - 



m + r to 



(16) 



(17) 



Effective potential is plotted in Figure^ for L = 3. The height of the 
potential decreases with the increase of the parameter u and finally merges 
with the Schwarzschild potential as u) tends to oo. 



3. Evaluation of quasinormal modes 

Third order WKB approximation method is used here to evaluate the 
complex normal mode frequencies of black hole, a semi analytic method 
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Figure 1: Effective potential for different values of to with L=3. Dotted line represents the 
Schwarzschild case. 



originally developed by Schutz and Willj32| for the lowest order. Later Iyer 
and Will [33] carried out this approach to third WKB order and Konoplya 34J 
to sixth order to get more accurate results. The method is found to be 
accurate for low lying modes and can be used to explore the QNM behavior of 
black holes efficiently, without going to the complicated numerical methods. 

WKB method gives a simple condition which will lead to discrete, complex 
values for the normal mode frequencies: 



k 2 = [V + (-2Vo')lA] - i(n + |)(-2T^)i(l + n), 
where A and Q are second and third order WKB correction terms, 



(18) 



A 



:-2^ ")v 2 




+ a' 



V 
v o 

288 V Vn 



(7 + 60a 2 ) 



(19) 



n 



1 

(-2Vo"; 
i 


.{ 5 

L 6912 


' 2304 


{< 


' 288 \ 





f) 4 ( 77+ 188a2 )-3^ 



(51 + 100a 2 ) 



(67 + 68a 2 ) + 

v ; 288 



(5 + 4a 2 )}, 



v"'v (5) 

y"2 



(19 + 28a 2 ) 



(20) 
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a 



n + | and n is the mode number, 



)( , 0,1,2, Re(E)yO 

1 -1,-2,-3, Im{E) -< 0. 1 > 

Vg denotes the n th derivative of V evaluated at r , the value of r at 
which V attains maximum. Tabled and Tabl42] list the QNMs evaluated. 





L-- 


=2 


L-- 


=3 


L- 


=4 


L-- 


=5 


0.5 


0.44322 - 


0.05945 i 


0.68305 - 


0.06558 i 


0.91050 


- 0.06764 i 


1.13252 - 


0.06863 i 


1 


0.40001 - 


0.07875 i 


0.63226 - 


0.08339 i 


0.84942 - 


0.08506 i 


1.06027 - 


0.08586 i 


1.5 


0.38949 - 


0.08221 i 


0.61995 - 


0.08694 i 


0.83446 - 


0.08857 i 


1.04250 - 


0.08933 i 


2 


0.38488 - 


0.08397 i 


0.61434 - 


0.08853 i 


0.82762 - 


0.09011 i 


1.03435 - 


0.09087 i 


.3 


0.38065 - 


0.08573 i 


0.60904 - 


0.09001 i 


0.82113 - 


0.09156 i 


1.02662 - 


0.09229 i 


4 


0.37866 - 


0.08661 i 


0.60650 - 


0.09072 i 


0.81801 - 


0.09224 i 


1.02290 - 


0.09297 i 


6 


0.37676 - 


0.08749 i 


0.60403 - 


0.09140 i 


0.81497 - 


0.09290 i 


1.01927 - 


0.09362 i 


8 


0.37583 - 


0.08792 i 


0.60281 - 


0.09174 i 


0.81348 - 


0.09323 i 


1.01748 - 


0.09394 i 


10 


0.37528 - 


0.08819 i 


0.60209 - 


0.09194 i 


0.81259 - 


0.09342 i 


1.01642 - 


0.09413 i 


oo 


0.37316 - 


0.08922 i 


0.59927 - 


0.09273 i 


0.80910 - 


0.09417 i 


1.02568 - 


0.09067 i 



Table 1: Fundamental^! = 0) QNMs of gravitational perturbations for various values of 
to and L. lj = oo represents the Schwarzschild limit 



n 


Ul- 


=1 


U)- 


=2 


IV- 


=3 


Ul- 


=4 




= oo 





0.63226 - 


0.08339 i 


0.61434 


- 0.08853 i 


0.60904 


- 0.09001 i 


0.60650 


- 0.09072 i 


0.59927 


- 0.09273 i 


1 


0.62029 - 


0.25228 i 


0.59996 - 


0.26827 i 


0.59381 - 


0.27290 i 


0.59085 - 


0.27511 i 


0.58235 


- 0.28141 i 


2 


0.59924 - 


0.42583 i 


0.57504 - 


0.45373 i 


0.56747 - 


0.46182 i 


0.56380 - 


0.46568 i 


0.55320 


- 0.47668 i 


3 


0.57207 - 


0.60378 i 


0.54310 - 


0.64417 i 


0.53372 - 


0.65589 i 


0.52913 - 


0.66149 i 


0.51575 


- 0.67743 i 


4 


0.54009 - 


0.78467 i 


0.50538 - 


0.83773 i 


0.49372 - 


0.85315 i 


0.48797 - 


0.86053 i 


0.47107 


- 0.88154 i 



Table 2: Higher mode frequencies for different values of u with L=3 

The variation of real and imaginary parts of QNMs with the parameter u 
are plotted in Figure^! The oscillation frequency, Re(u) decreases with the 
increase of w while the damping time \Im(uj)\, increases with the increase 
of w. For values near w = 0.5 there is a sudden decrease in damping time. 
Finally comparing with the Schwarzschild black hole, oscillation frequencies 
of QNMs in HL gravity is larger and has a lower damping time. 

4. Conclusion 

We have investigated the evolution of gravitational perturbations in black 
hole spacetime in HL theory and associated QNMs are evaluated. The nega- 
tive imaginary part of QNMs show that the perturbation will decay in time 
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Figure 2: Re(k) and Im(k) as a function of w for n = and L=2,3,4,5. 

and hence we can conclude that black holes are stable against small perturba- 
tions in the spacetime. The present calculations show that the gravitational 
QNMs are long lived in Hof ava theory compared with the Schwarzschild black 
hole. 
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